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Abstract. Mode-coupling theory for uniformly sheared underdamped systems with an isothermal condition is presented. As 
a result of the isothermal condition, it is shown that the shear stress exhibits significant relaxation at the a-relaxation regime 
due to the cooling effect which is accompanied by the growth of the current fluctuation. This indicates that nonequilibrium 
underdamped MCT is not equivalent to the corresponding overdamped MCT, even at long time-scales after the early stage of 
the jS-relaxation. It is verified that the effect of correlations to the density-current modes is, however, negligible in sheared 
thermostated underdamped systems. 

Keywords: glass, mode-coupling theory, underdamped systems, response 
PACS: 64.70.P-, 61.20.Lc, 83.50.Ax, 83.60.Fg 



INTRODUCTION 

Glassy liquids are extensively studied for understand- 
ing the transition from liquid to amorphous solid states. 
Two typical theoretical approaches coexist in these stud- 
ies; the mode-coupling theory (MCT) [1], and the replica 
theory with assumptions on the metastable states from 
the random-first-order transition theory (RFOT) [2], 
Though not yet complete, the relation between the above 
two approaches has been revealed in some aspects [2, 3]. 
For instance, it has been shown that the equation for the 
non-ergodic parameter (NEP) in MCT can be partly de- 
rived by the replica theory [2] . However, since the replica 
theory is a static theory, it lacks genuinely dynamical ef- 
fects, and its application to the study of dynamics, e.g. 
relaxation towards steady state, or response, is difficult. 

MCT, which is a dynamical theory, also suffers from 
difficulties when applied to glassy liquids. For instance, 
it predicts unphysical non-ergodic transition, which can- 
not be observed in realistic thermal systems [4, 5]. 
Another issue is the lack of a consistent perturbative 
field-theoretic formulation, which is practically applica- 
ble to higher orders. Although efforts have been made 
[6, 7, 8, 9], they are still not comparable to the con- 
ventional approach of projection operators and the Mori- 
type equations [10], which we adopt in this work. 

Glassy liquids can be classified into two categories; 
one is molecular glasses, which is described by New- 
tonian dynamics, and the other is colloidal glasses, 
which is described by Brownian dynamics, at the micro- 
scopic level, respectively. Historically, MCT for molec- 
ular glasses has been developed, and has passed various 
tests concerning the two-step relaxation phenomena [1]. 
Then, it was shown theoretically that the MCT for col- 



loidal glasses are equivalent to that for the molecular 
glasses in the long-time limit [11]. That is, the height 
of the plateau of the density time-correlator at the /3- 
relaxation regime (NEP), and the a-relaxation time, are 
coincident. This fact has been verified by molecular dy- 
namics (MD) simulations [12]. Since then, it has been 
believed that the glassy behavior is insensitive to its mi- 
croscopic origin. 

However, the established equivalency between the two 
categories is only verified for equilibrium MCT, and is 
not yet established for nonequilibrium MCT. It is non- 
trivial whether nonequilibrium MCT for underdamped 
systems, such as sheared molecular liquids or granular 
matters, share equivalent long-time behaviors with that 
for overdamped systems. 

Nonequilibrium MCT for uniformly sheared over- 
damped systems, for e.g. colloidal suspensions immersed 
in solvents, has been worked out by the pioneering pa- 
pers of Fuchs and Cates [13, 14] and Miyazaki et al. 
[15, 16]. In these works, it has been shown that the shear 
thinning and the emergence of the yield stress in glassy 
states can be captured within this framework. Further- 
more, comparison of MCT with MD has been performed 
in Ref. [16], and it has been shown that, aside from 
the difference in the early stage before the /3 -relaxation 
due to inertia effects, the plateau of the density time- 
correlator coincide with each other. 

On the other hand, nonequilibrium MCT for uniformly 
sheared underdamped systems has been formulated by 
Chong and Kim [17], whose microscopic dynamics is 
described by the SLLOD equations [18]. However, there 
were some problems in this formulation, as indicated in 
Ref. [19], and calculations to verify its validity have not 
yet been performed. In contrast to equilibrium MCT and 



sheared overdamped MCT, where the temperature is a 
mere parameter, the temperature is a dynamic variable 
in sheared underdamped MCT. Hence, it is necessary to 
control the temperature to avoid heat-up or cool-down in 
order to conform with realistic physical situations. One 
of the problems of Ref. [17] is that the temperature is not 
controlled and heats up. 

In this paper, we formulate the MCT for uniformly 
sheared underdamped systems with an isothermal con- 
dition imposed. We present that there is a significant re- 
laxation in the shear stress at the a-relaxation regime, 
although the profile of the density time-correlator is coin- 
cident with the overdamped systems after the early stage 
of the j3 -relaxation. Then, we address another issue spe- 
cific to underdamped systems: the significance of corre- 
lations to the density-current modes. We show that they 
are negligible in sheared underdamped thermostated sys- 
tems, but discuss that this situation might be acciden- 
tal. This result is new, which has not been discussed in 
Ref. [19]. 

We stress that this paper is related to our previous pa- 
per [19] in some respects, but it addresses issues not in- 
cluded there, such as the effect of the density-current 
correlation, as well. We also point out that our formu- 
lation of the underdamped MCT developed in this paper 
can resolve the problem of the stress overshoot in the 
a-relaxation regime [20, 21]. It is shown that the over- 
damped MCT underestimates the stress overshoot com- 
pared to MD [20], while it is implied that it is consistent 
with the Brownian dynamics simulations [21]. Thus, it 
might well be concluded that the pronounced stress over- 
shoot in molecular liquids than in colloidal suspensions 
is the result of the inertia effect. 



UNDERDAMPED ISOTHERMAL MCT 

In this section, we present the MCT of uniformly 
sheared underdamped systems, with isothermal condi- 
tion imposed. The basis of the formulation is the mi- 
croscopic SLLOD equations [18] for a spherical N-body 
system of mass m, 

Ht) = Z^ + x-rtit), (1) 
m 

p,(t) = Ft(t)-K- Pi (t)-a(r)p t {t), (2) 

where i = l,--,N, and their corresponding Liouville 
equations. Here, K Xfl — fS^S^ = x,y,z) is the 

shear-rate tensor, where y is the constant shear rate, and 
Fj(t) is a conservative force with a soft-core potential 
whose details are irrelevant in the following discussions. 
The parameter a(T) in Eq. (2), where T = {ri 7 Pj}f =l 
with r = r(0), gives the dissipative coupling of the par- 
ticle with the thermostat. If we impose the isothermal 



condition, which constrains the kinetic temperature to re- 
main constant from its initial equilibrium value, <x(V) 
should be regarded as a multiplier for this constraint, 
rather than an independent physical parameter. A repre- 
sentative example of this is the Gaussian isokinetic (GIK) 
thermostat [22], which is established as one of the most 
typical thermostat in MD. The specific form of a(T) in 
this thermostat is 

g(r) = a,(^„-WrO 

which follows from the constraint 4- Y^ =l Pi(t) 2 = 0. 

Before the discussion of the isothermal condition 
in MCT, we address the basics of the sheared under- 
damped MCT. The major difference of the underdamped 
to the overdamped MCT is that, besides the density 
time-correlator ^> q (t), we should consider the current- 
density cross time-correlator H q (t) (X = x,y,z) on the 
same grounds. A crucial feature of the sheared system is 
that translational invariance is preserved in the sheared 
frame, although it is broken in the experimental frame. 
Hence, Fourier transform in the sheared frame is valid. 
This leads to the following selection rule of wavevec- 
tors for two-point functions, (T(f ))B q (T)*')^ = 

(A q(t) (r(t))B q (r)*)^8 k , q . Here, q(t) = q - (q ■ K)t is 
the Affine-deformed wavevector, which incorporates the 
effect of shearing. This selection rule leads to the follow- 
ing definitions of the above two time-correlators, 

**(0 = ^(" ? w( f K(o)*) eq , (4) 

H tt) s ^4w( f K(or) eq , ( 5 ) 

where n q {t) = - N8 qfi and 4(f) = 

Yd=i pf (t)e' q ' r '^ jm are the density and the current- 
density fluctuations, respectively. These definitions are 
not equivalent to those of Ref. [17], where Fourier 
transform is performed in the experimental frame. 

Now we discuss the implementation of the isothermal 
condition. The most straightforward and rigorous way to 
perform this is to adopt as a (F) the one for the GIK ther- 
mostat, Eq. (3). If this is possible, the resulting Liouville 
and Mori-type equations are expected to automatically 
satisfy the isothermal condition. To rephrase it, not only 
the averaged temperature, but also the microscopic tem- 
perature, is expected to remain constant. However, this 
does not fit MCT, because the integral of a rational func- 
tion of momentum, Eq. (3), with Gaussian weight is dif- 
ficult to perform explicitly. Hence, we propose another 
way, which only requires the averaged temperature to re- 
main constant. This is attained by requiring the isother- 
mal condition to be satisfied at the level of the Mori-type 
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FIGURE 1. Time-evolution of the normalized density time-correlator <fr q (t)/S q . The left (right) figure is for the underdamped 
(overdamped) systems. Time is scaled in unit of To = d/vo, where d is the diameter of the spherical particle and vo is the relative 
shear velocity at the boundaries. The three lines correspond to nondimensionalized wavenumbers qd = 5.0, 6.2, 7.0. The shear rate 
7 is chosen as yTo = 10~ 4 , and the volume fraction <p is chosen as e = (<p — <p c )/<Pc — +10~ 3 , where (p c ~ 0.516 is the critical 
volume fraction of the equilibrium MCT transition. 



equation, which is an equation for the averaged time- 
correlator. 

The averaged temperature is defined by (r(f)) eq = 
2{K(r(t))) eq /(3Nk B ), where K(T) ee lliPf/(2m) is 
the total kinetic energy, and (■ ■ •) represents the equi- 
librium ensemble average. The time derivative of the av- 
eraged temperature can be obtained from the generalized 
Green-Kubo formula [23], and the isothermal condition 
reads 



d , . 2 



(^(r(0)a(T(o))) Bq = o ) (6) 



where 



n(T) = -j8V7^(r)-2j8a(r)5i:(r) 



(7) 



is the work function. Here, V is the volume of the system, 
j3 = l/(ksT e q) is the inverse equilibrium temperature, 
8K(T) is the fluctuation of the kinetic energy, and a(T) 
is a multiplier for the constraint Eq. (6) which appears 
in the Liouvillian, whose explicit form is discussed be- 
low. The work function £l(T) represents the gross work 
loaded on the system by shearing and the dissipative cou- 
pling to the thermostat. 

By introducing general time-correlators of the form 

G A , B (t) = <A(r(0)B(r)) eq = ([u(t)A(r)}B(r)) ea , (8) 



cq : 



where U(t) is the time-evolution operator, the isother- 
mal condition Eq. (6) is written in terms of these time- 
correlators as 



YG K> a(t)+2G K „ SK (t) = 0. 



(9) 



This condition indicates a balance between the aver- 
age work by shearing and the average dissipation. We 
introduce a multiplier for the constraint Eq. (9) in 
G K a g K (t), which we denote by X a (t). This is attained 
by defining a "renormalized" time-evolution operator 
J/#(f) = k a {t)U{t), and a corresponding time-correlator 



T K,aSK 



(f), where 



([U R {t)K{T)]a{T)8K{T)) eq 
K{t)G K „ 5K {t\ (10) 



with which we recast Eq. (9) as 

YG K Mt)+2G^ aSK (t)=0. 



(11) 



Note that X a (t) is a multiplier which is averaged and ap- 
pears in the Mori-type equations, and hence independent 
of r. From Eqs. (10) and (1 1), we can see that the isother- 
mal condition is satisfied by choosing X a (t) as 



X a (t) 



7 G K .a{t) 

2 G KM g K (t) 



(12) 



if G K ^ a $ K (t) ^ 0. On the other hand, if G K a sK(t) = 0, 
the system lacks cooling, and it inevitably heats up. 

The Mori-type equations for the correlators <&q(t) and 
Hy (t ) are derived by introducing the projection operator, 



E 



Xn 



eq 



NS, 



n k( t ) - 



Xj k(t), 



k Nv 2 T 



Here, S# is the static structure factor and vj = \jT &q /m is 
the thermal velocity, where r eq is the initial equilibrium 




FIGURE 2. Time-evolution of the average shear stress (left) and the longitudinal component of the current-density cross time- 
correlator Hq(t) (right). In the underdamped case, the shear stress is significantly relaxed at the a-relaxation regime, which is 

accompanied by a growth of H^(t). The shear rate is chosen as j^o = 10~ 4 , and the volume fraction as £ = (<p — <p c )/<J) c = +10~ 3 . 



kinetic temperature. To proceed further, it is necessary 
to resort to the mode-coupling approximation, which 
determines the specific choice of 05(1") and X a (t), and 
provides a closure of the Mori-type equations. To this 
end, the second projecton operator to the pair-density 
modes 

Y^ n *k{t) n *p(t 

^nn{t)X = £ -^—i^ n m n p(f) (14) 

is introduced. The validity of this choice will be dis- 
cussed in the next section. The simplest form of (x(T) 
is a(T) = Oo = const., where no fluctuation is incorpo- 
rated. This is the choice adopted in Ref. [17]. However, 
we have figured out that this choice leads to G K a g K (t) = 
0; the isothermal condition cannot be satisfied, and the 
system heats up. Hence, it is necessary to incorporate 
fluctuations into ce(T). In this respect, the simplest form 
would be 



«(T) 



^Nk B T eq ^2m 7 



(15) 



where current fluctuations are incorporated. Here, Oo 
is a constant. With this choice, it can be verified that 
GK.aSK^) 7^ 0' ar, d hence the constraint Eq. (11) can be 
satisfied. Straightforward calculation of G K a § K (t) and 
Gk aif) leads to 



A a (f)ao 



y£*>° s m s k 



W k <S> k (tf 



2 Lie 



>0 s m s k 



<S> k (t) 2 



(16) 



where W k = ^x~3^1?F t^]- This expression corresponds 
to a(T) of the GIK thermostat; it is a ratio of the work 



loaded on the system (numerator) to the total kinetic 
energy (denominator). 

Now we present the Mori-type equation with the 
mode-coupling approximation applied. We adopt the 
isotropic approximation, in which case the two Mori- 
type equations can be cast in a single second-order equa- 
tion for & q (t), 



dt 



l a (t)ao-f- 



m 2 



s, 



%(t)- I dsM q{s) (t-s)-%(s), (17) 



where the multiplier X a (t) appears as an effective fric- 
tion coefficient. It is notable that the memory kernel 
is unaltered by the introduction of the multiplier with 
the choice of the second projection operator Eq. (14). 
In the isotropic approximation, only the longitudinal 
component of the current-density cross time-correlator, 
Hq(t) = d<P c/ (t)/dt, is considered. The numerical solu- 
tion of Eq. (17) is shown in Fig. 1, together with the 
corresponding result for the overdamped case. As can be 
seen, the long-time behavior after the early stage of the 
j3 -relaxation, namely the height of the plateau of the den- 
sity time-correlator and the a-relaxation time, is equiva- 
lent for the two cases. 

Despite this fact, there is a significant difference be- 
tween underdamped and overdamped MCT for sheared 
systems. From the generalized Green-Kubo formula, the 
following formula for the shear stress can be derived for 
the underdamped case [19], 



<<*(')> 



k R T t 



cq 



ds 



d 3 k 



(18) 



The corresponding well-known formula for the over- 
damped case [13] can be obtained by setting X a (t)(Xo =0 
in Eq. (18). This term with X a (t)(Xo is exactly the "cool- 
ing effect" specific for sheared underdamped systems, 
which prevents the system from heating up. The numer- 
ical result of Eq. (18), together with that of H^(t), is ex- 
hibited in Fig. 2. We can see that, in the underdamped 
case, there is a significant relaxation of the shear stress 
in the a-relaxation regime. Furthermore, this relaxation 
is accompanied by a growth of H^(t). From these re- 
sults, we can obtain an intuitive picture as follows; during 
the j3 -relaxation regime, the particles are caged by their 
neighbors, and they are almost frozen. When the accu- 
mulated strain by shearing becomes comparable to the 
particle size, the cage is broken, and the partiles start to 
escape from the cage. This generates current fluctuations, 
which is cooled to keep the temperature of the system un- 
changed. This cooling effect leads to a relaxation of the 
system to the genuine stable state, and a stress relaxation 
occurs. 

It is remarkable that the overdamped MCT underes- 
timates the stress overshoot in the a-relaxation regime 
compared to its corresponding MD [20]. On the other 
hand, it is impled that the overdamped MCT predicts the 
magnitude of the stress overshoot consistently with the 
Brownian dynamics simulations [21]. This consistency 
is established by means of the schematic model, rather 
than the microscopic model, of the MCT, so a complete 
verification remains to be a future problem. However, we 
might say that the results of Refs. [20] and [21] strongly 
suggest that the pronounced stress overshoot in molecu- 
lar liquids than in colloidal suspensions is the result of 
the inertia effect. 



EFFECT OF DENSITY-CURRENT 
MODES 

In the previous section, we have addressed a prob- 
lem specific to the sheared underamped MCT; that is, 
to control the temperature from heating-up or cooling- 
down. There is another non-trivial problem specific to 
the sheared underdamped MCT which should be ad- 
dressed; "up to what degree should we take correla- 
tions into account?". The fundamental concept of MCT 
is to take into account correlations with "slow" modes. 
In the overdamped case, it is the density fluctuation 
n q {t), and in the underdamped case, it is n q (t) and 
the current-density fluctuation j q (t). Although jg(t) is 
"faster" than n q (t), it is necessary to incorporate to deal 
with physical phenomena such as stress relaxation at the 
a-relaxation regime, which we have shown in the pre- 



vious section. Incorporating correlations, or projecting 
onto n q (t) and/or jg (f ), leads to the Mori-type equations, 
which includes a memory kernel as a time-correlator of 
noises (uncorrected degrees of freedom). A memory ker- 
nel includes higher-order correlations, and we should ex- 
tract these correlations by projecting onto higher-order 
modes, such as nn, n j, nnn, (Of course, correlations 
with higher-order modes correspond to multi-body cor- 
relations, which must be truncated to close the equation. 
MCT is a scheme which approximates a four-body cor- 
relation by factorizing it into a product of two-body cor- 
relations.) For this purpose, we have introduced the sec- 
ond projection operator Eq. (14), which is conventional 
in equilibrium as well as in sheared overdamped MCT. 
However, in the sheared underdamped MCT, it is not 
a priori apparent that the projection onto the density- 
current modes 



E 

k, P 



Xn 



k(t)jp(t 



N 2 S, 



k(t) v T 



n mJp(t) 



(19) 



is negligible compared to £? nn (/), and hence the adoption 
of Eq. (14) must be validated. In particular, we should 
carefully check the roles of Z^nji?) as well as the self- 
consistency of our argument in the previous section in 
our underamped MCT, because the existence of the mo- 
mentum current causes some essential differences in the 
long-time dynamics from that in the overdamped MCT. 
In this section, we discuss the effect of the projection 
onto the density-current modes in the second projection 
operator, Eq. (20). 

Let us consider a second projection operator 2?i{t) of 
the following form [24], 



(20) 



where 8P nn (t) and & n j(t) are defined by Eqs.(14) and 
(20), respectively. We note that the isothermal condi- 
tion Eq. (11) and the stress formula Eq. (18) are un- 
altered by the introduction of £P n j(t). Crucial obser- 
vation for this property is that the following equal- 
ities hold, ^njl^OxyiT) = 0, &> n] (t)K{r) = 0, and 
3» n j(t) [a{r)SK(r)] = 0, which follow from the vanish- 
ing of ensemble averages with odd number of momen- 
tum variables. Here, the explicit form Eq. (15) of a(T) 
is assumed. As a result, the time-correlators Gic,a{t), 
G K,aSK(t), G a , a (t), and G a , aSK (t) remain unchanged. 
The isothermal condition and the stress formula are ex- 
pressed in terms of G K aSK (t), G a aSK (t) and G K ,a{t), 
Ga,a(t), respectively, so the above statement is proved. 

Corrections due to S? n j{f) arise in the Mori-type equa- 
tions for the time-correlators. The introduction of 2? n j{t) 
requires us to introduce two additional time-correlators 
defined as follows, 



N 



(21) 




FIGURE 3. Numerical results for the time-correlators 4> 9 (f) (left) and 'Pg(t) (right). In the left figure, solid and dashed lines 
with caption (nj) are the solutions of Eq. (25) with ^2(0 = ^nn(t) + &nj(t), while circles and crosses with caption (nn) are the 
solutions for @>i{t) = 3? n n(t). 



(22) 



in addition to the previously introduced (t) and (f ), 
which are defined in Eqs. (4) and (5), respectively. Note 
that Hq (?) and H q (t) are inequivalent, since sheared sys- 
tems do not possess translational and reversal symme- 
tries in time. Similarly to ^q(t) and H^(t), we adopt the 

isotropic approximation for (t ) and (t). By intro- 
ducing two scalar time-correlators, ^ q {t) and C q (t), it 
reads 



-q(t) l V q (t), 



C*"(0 ~ S x *C q (t). 



(23) 
(24) 



By this approximation, it is possible to derive the follow- 
ing equations for the time-correlators, 



-[ka(t)ao-2Mlr,t)]-<t> q (t) 
-J o dsM q{s) {t-s)—%(s), (25) 



dt 2 



V 2 T q{tf 



+ A a (f)oo/i (j,t) 



3 S q{l) 

[A a (f)ao + /i(7,0]^%(0 



3 Jo 



, (26) 



where functions fi(f,t) and f2{j,t) which incorporate 
anisotropic effects of shearing are defined as 



MM 
MM 



1 <2- f 2\ 



(1 + I7V) 



(27) 
(28) 



Here, we have eliminated C q (t) by the relation 

C q (t) = 



d , N 1 c/ , n7 



(29) 



which follows from the Mori-type equation. The memory 
kernel in Eq. (25) consists of two terms, which are given 
by 



M q (x) = £ M«(t), 

i=l,2 

,,(1)/ v «Vj r 

* = 5# 



(2ff)3 



(9 • *) c Mt) + (9 • P)cp(T) [(q ■ k)c k + (q ■ p)c p ] 



x4> i: (T)O p (T), 



(30) 



M ( q 2 \x) 



2A„(t)0£o 



1 , ^' 

1 + 3 (P) 



d 3 k 



{<l-k)c-u T) + {q-p)c p(z) (q p 



x (t) 



(31) 



Note that the first term M q ' (t) is nothing but the conven- 
tional memoery kernel familiar in equilibrium as well as 
in sheared overdamped MCT, which is the source of the 



plateau in the density time-correlator. The second term 

(2) 

Mq (t) is the novel therm which originates in the "cool- 
ing effect" due to the isothermal condition. The memory 

kernel in Eq. (26), which is the trace of Mq^(x), also 
consists of two terms, which are given by 



(i)ll 



M, 



(1)11 



i=l,2 

nv\ r d 3 k 
x4> A; (t)4> p (t 



[kc k + pc p ] ■ 



M 



(2)11 (T) = 



-2A a (T)oco 
xp- 



1 



I 



d 3 k 



PCp(r 

(32) 



feq ( T)+p^ (T) l<f fc (T)^(T), (33) 



where again the second term originates in the "cool- 
ing effect". The two correlators <E> ? (?) and *¥ q {t) cou- 
ple through M { q 2) (x) and M ( q 2)XX (x). The derivation of 
Eqs. (25)-(33) is somewhat lengthy, though straightfor- 
ward, so we present only the results. 

It should be noted that, in the absense of the coupling 
to the thermostat, M^\z) and M^ XX (x) vanish. This 
suggests that the memory kernel is coincident for the 
equlibrium underdamped and overdamped MCT, even 
when projection onto the density-current modes ^ n j(t) 
is included. This proves the validation of the projection 
onto the pair-density modes 2? nn (t) in the equilibrium 
underdamped MCT. 

The numerical solutions of Eqs. (25) and (26) are 
shown in Fig. 3. For the density time-correlator $> q (t), 
the results for the case without including ZP n j(t), i.e. 
2?i{t) = &nn(t)> ^ plotted in circles and crosses. The 
initial conditions are chosen as follows: 



%(t = 0) = S q , 

«F,(f = 0) = 0, 



d 
dt 

5™ 



-0, (34) 



r=0 



r=0 



^. 05) 



As can be seen from Fig. 3, although a sign of corre- 
lation of currents can be seen in VP 9 (f) (which appears 
as a peak around ?/to — 10 _1 ), the result of <t> q (t) is 
almost coincident with that for 2?2{t) = £P nn (t). To- 
gether with the fact that ^(t) of Eq. (20) leads to the 
identical isothermal condition and the stress formula as 
for ^2(0 = ^nn{t), it is verified that projection onto 
the density-current modes is negligible in sheared ther- 
mostated systems as well. Thus, the analysis in the pre- 
vious section, where we only use £? nn {f) for the second 
projection operator, is self-consistent. 

In principle, we can see from Eqs. (31) and (33) 
that the "cooling effect" which affects the memory ker- 



nel works to destroy the plateau of the density time- 
correlator, since their contributions are negative. This can 
be interpreted as the exhausion of the cage by cooling, 
which allows the caged particles to escape. However, as 
shown in Fig. 2, Mq 2 \x) and Mq 2 ^ X (x) becomes signif- 
icant only at the a-relaxation regime, where the plateau- 
generating kernels M q l \x) and Mq l ^ X (x) decay irre- 
spective of the existence of Af^(r) and (t). This 

is why the effect of (t) and M^ 11 (x) is hidden be- 
hind. 

Although the role of 3? n j(t) is not significant in the 
present case, it is an important feature that dissipation 
acts to destroy the plateau of the density time-correlator. 
A striking example can be found in granular liquids [25], 
where the plateau of the density time-correlator is de- 
stroyed by inelastic collisions. Hence, we might rather 
consider that correlations with density-current modes 
should be included in general, and it might be accidental 
that they can be neglected in sheared thermostated sys- 
tems. 

Note that the result of *P 9 (f) (which is essentially 
Hq{t)) does not show a growth in the a-relaxation 
regime, in contrast to 2i^q{t) (which is essentially 
is an explicit manifestation of the inequiva- 
lence of h£ (t) and (f ) mentioned previously. 



DISCUSSION AND CONCLUSION 

At the formal level, response formulas can be derived 
from the generaized Green-Kubo (GGK) relation [23, 
26]. For the case of constant external fields, this relation 
gives a formula for the steady-state averages, which char- 
acterizes the relaxation towards a steady state. In con- 
trast, for the case of time-dependent external fields, the 
time-ordering of operators should be handled with care, 
whose application to the GGK relation is relatively com- 
plicated. As an alternative, we have constructed an ex- 
plicit approximate (quasi-)steady-state distribution func- 
tion valid on the plaetau of the density time-correlator, 
and have derived a response formula at the linear order. 
The situation we consider might be restricted, but at least 
it is explicitly calculable, and can be compared to exper- 
iments. 

In conclusion, we have formulated a nonequilibrium 
MCT for uniformly sheared underdamped systems with 
an isothermal condition. We have figured out that there is 
a significant stress relaxation at the a-relaxation regime 
due to a cooling effect, which is accompanied by the 
growth of current fluctuations. This stress relaxation is 
predicted to be larger than that in the overdamped MCT, 
which is consistent with MD and Brownian dynamics 
simulations [20, 21]. This indicates that nonequilibrium 



underdamped MCT is not equivalent to its correspond- 
ing overdamped MCT, even at long time-scales after the 
early stage of the j3 -relaxation. This is in sharp contrast 
to the equilibrium case. We have also shown that corre- 
lations with density-current modes are not important in 
sheared underdamped thermostated systems, but argued 
that this might be accidental, and they should be incoro- 
prated in general. 

Another example which cannot be described by the 
overdamped formulation is granular particles [25]. The 
study of rheological properties of granular liquids by 
means of sheared underdamped MCT presented here is 
in progress, and will be reported elsewhere [24]. 
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